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1. Introduction 

Generalized parton distributions (GPDs), introduced some time ago, 1-3 
have received much attention from both the theoretical and the experimen- 
tal side. This was triggered by the 'spin crisis', referring to the mismatch of 
quark spin, extracted from polarized deep inelastic scattering, and the con- 
stituent quark model. This is rather analogous to the 'momentum crisis', 
where almost 50% of the proton momentum is carried by gluons. Obviously, 
both 'crises' of the constituent quark model are 'solved' by switching to the 
fundamental degrees of freedom. Hence, we view the 'spin puzzle' first as 
a quest to quantify the partonic structure of the nucleon in terms of quark 
and gluon angular momenta and second as a wish for an appropriate match- 



*Talk given by D.M. at the 4th Workshop on Exclusive Reactions at High Momentum 
Transfer, May 18-21, 2010, Thomas Jefferson National Accelerator Facility, Newport 
News, Virginia USA. 
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ing of effective and partonic degrees of freedom. Furthermore, it has been 
realized that GPDs build up a whole framework for description of hadron 
structure, 4 ' 5 including a three dimensional imaging. 6,7 Certainly, the 'spin 
puzzle' is one interesting and intricate aspect, which we will now discuss. 

An appropriate decomposition of the nucleon spin in terms of energy- 
momentum tensor form factors has been suggested by X. Ji: 8 

l - = jQ + J G , JQ= £ J\ r = A*(Q 2 ) + B*(Q 2 ), (1) 

q—u,d,--- 

where the quark and gluon contributions are given by the first moments of 
parity even and target helicity (non-)conserved GPD H (E), 

{ ^ }(Q 2 ) = Jjm^ fdx (x, r,, t, Q 2 ) , i E {u, d, ■ ■ ■ , G} ,(2) 

taken in the forward limit. Note that the momentum fractions A 1 are ob- 
tained from global PDF fits and that momentum and angular momentum 
conservation imply that the total gravitomagnetic moment vanishes, 

J2 B«(Q 2 ) + B G (Q 2 )=0. 

q—u,d^--- 

Hence, the unknowns in the spin sum rule are the helicity flip form factors 
B l . Any 'pure' quark model trivially predicts A G = and B G = 0, whereas 
taking for A G the phenomenological value ~ 0.45 and relying on B G = 0, 
one concludes that quarks and gluons carry roughly the same amount of 
angular momentum a . Indeed, lattice calculations indicate a rather small 
B®. Unfortunately, these simulations suffer from systematic uncertainties, 
in particular, so-called disconnected contributions are presently neglected. 10 
Thus, a phenomenological handle on the spin sum rule is highly desired. 

Deeply virtual Compton scattering (DVCS) off nucleon is considered as 
the theoretically cleanest process offering access to GPDs. Its amplitude 
can be parameterized by twelve Compton form factors (CFFs), 11 which are 
given in terms of GPDs, e.g., at leading order (LO) in generic notation: 

{?}(«. t. o") L =° /> { I }(*,,- f , t, a 2 ) • (3) 

The Bjorken variable xb might be set equal to 2£/(l + £). Favorably, DVCS 
enters as a subprocess into the hard photon electroproduction where its in- 
terference with the Bethe-Heitler bremsstrahlung process provides variety 



a The orbital angular momentum of quarks is defined as L q = J q — S 9 /2, where E 9 /2 ~ 
0.15 is the spin of quarks. For more detailed discussions see review. 9 
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of handles on CFFs. However, the target helicity conserved CFF H is of- 
ten the dominant contribution, while £ always appears with a kinematic 
suppression factor t/AMfj, induced by the helicity flip. In other words, one 
should look at observables for which % is suppressed, too, which requires a 
transversely polarized target. 12 Also, in photon electroproduction off neu- 
tron 13 Ti is suppressed in the interference term by the accompanying Dirac 
form factor F™ (F™(t = 0) =0). Unfortunately, one also has to worry about 
other non-dominant CFF contributions. Thus, the extraction of £ requires 
a most complete measurement of all possible observables in dedicated ex- 
periments. 

Supposing £ is measured, the question arises: Can one deconvolute 
Eq. (3)? Apart from radiative and higher twist-contributions, one might 
view the GPD on the r\ = x cross-over line as a "spectral function" , which 
provides also the real part of the CFF via a dispersion relation: 14-17 

3mJ( lB ,i,S 2 ) L =° 7rF(£,£,i,Q 2 ), F = {H,E,H,E} 7 (4) 
3fe{^}(su,t, Q 2 )^Pv£dx^^{ H E y x ,x,t,Q 2 )±V(t, Q 2 ).(5) 

The GPD support properties ensure that Eqs. (4,5) are in one-to-one corre- 
spondence to the perturbative formula (3), where the subtraction constant 
T> can be calculated from either H or E. To pin down the GPD in the outer 
region y > i] = x, one might employ evolution, e.g., in the non-singlet case, 

V 2 -^2F(x,x,t,V 2 )= I ^V(l,x/y,a s (^)F{y 7 x 7 t^ 2 ). (6) 
a/j, J x x 

A large enough Q 2 range is not available in fixed target experiments. Hence, 
we must conclude that essentially only the GPD on the cross-over line 
(thanks to (5), also outside of the experimentally accessible part of this 
line, 17 ) and the so-called £>-term 18 can be accessed. Moments, such as those 
entering the spin sum rule (1), can only be obtained from a GPD model, 
fitted to data, or more generally with help of some 'holographic' mapping: 17 

{F(x,r ] = 0,t,Q 2 ),F(x,x,t,Q 2 )} =► F(x,r],t, Q 2 ) . (7) 

Here, F l (x,n = 0, i, Q 2 ) are constrained from form factor measurements 
and, additionally, GPDs H % (W) by (un)polarized phenomenological PDFs. 
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2. GPD representations and modeling 

Let us now turn to GPD modeling in different representations. First, GPDs 
might be defined as Radon transform of double distributions: 1,19 (DD) 



where p 6 {0, 1}. In this representation polynomiality, however, not pos- 
itivity constraints 20 are explicitly implemented. Moreover, with the right 
choice for p, h the polynomiality of ir-momcnts can be completed to the re- 
quired order in r\. In the central, —r\ < x < i], and outer, r\ < x < 1, region 
the GPD can be interpreted as the probability amplitude of a i-channel 
meson-like and s-channcl parton exchange, respectively. Mathematically, F 
is a twofold image of the DD /, where those in the central and the outer 
region can be mapped to each other. 15,17 ' 21 The potential ambiguity, a term 
that lives only in the central region, is removed by requiring analyticity. 15,17 
Popular GPD models are based on Radyushkin's DD ansatz 19 for t = 0, 
where the DD factorizes into the PDF analogue f(y) and a normalized 
profile function II(z). The GPD on the cross-over line is then given as 



which is a linear integral equation of the first kind within the kernel 
/( ^'^ /(l — xz). Knowing the GPD at r\ = 0, i.e., f{y), and on the 
cross-over line, allows to determine the profile function and so to recon- 
struct the entire GPD C , giving example of the 'holographic' mapping (7). 

On the first glance a GPD in the outer region can be straightforwardly 
represented by an overlap of light-cone wave functions (LCWF), 23 which 
guarantees that positivity constraints are implemented. In simple models 
one might even reduce the number of non-perturbative functions, e.g., in a 
spectator diquark model, one only deals with one effective scalar LCWF for 
each struck quark species. This predicts for each of them four chiral even 
and four chiral odd GPDs. Moreover, one might use such a representa- 
tion to evaluate also transverse momentum dependent parton distributions 
(TMDs). Unfortunately, there is a drawback. In the central region the GPD 
possesses an overlap representation in which the parton number is not con- 
served, and where the LCWFs are dynamically tied to those used in the 



b Note that the factor (1 — x) might be replaced by a more general first order polynomial. 
c An example is provided by f(x) oc x~ a (l — x)@ , which yields after some redefinitions 
the integral kernel k = (1 — rrz) - ' 9-1 . The solution is then obtained in Mellin space. 22 
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outer region. A closer look reveals that Lorentz covariance already ties the 
momentum fraction and transverse momentum dependence of a LCWF, 24 
see also Refs. 25-27. Hence, a overlap representation is only usable if the 
LCWFs respect Lorentz symmetry, which would allow to restore the GPD 
in the central region. 24 

Strictly spoken, positivity constraints for GPDs are only valid at LO, 
since they can be violated by the factorization scheme ambiguity. Never- 
theless, it would be desired to impose them on GPD models. One might 
follow the suggestion 28 and model GPDs as an integral transform of (tri- 
angle) Feynman diagrams, i.e., spectator quark models. A specific integral 
transformation, namely, a convolution with a spectator mass spectral func- 
tion, can be used to include Regge behavior from the s-channel view. Such 
dynamical models provide also effective LCWFs or TMDs; however, sim- 
plicity is lost. In particular, PDF and form factor constraints cannot be 
implemented, i.e., one has to pin down such models within global fitting. 

At present we neglect positivity constraints and model GPDs in the 
most convenient manner by means of a conformal SL(2,R) partial wave 
expansion, which might be written as a Mellin-Barnes integral 21 

Here, pj(x,r]) are the partial waves, given in terms of associated Legendre 
functions of the first and second kind, and the integral conformal GPD 
moments Fj(r),t, /z 2 ) are even polynomials in rj of order j or j + 1. The 
advantages of this representation are: i. the conformal moments evolve au- 
tonomously at LO, ii. one can employ conformal symmetry to obtain NNLO 
corrections to the DVCS amplitude 15,29 and iii. PDF and form factor con- 
straints can be straightforwardly implemented. Namely, Fj (r] = 0, t = 0, tt 2 ) 
are the Mellin moments of PDFs, Fj =0 are partonic contributions to clastic 
form factors, Hj = \ and Ej = \ are the energy-momentum tensor form factors, 
and for general j one immediately makes contact to lattice measurements. 
One might expand the conformal moments in terms of t-channel SO(3) 
partial waves 30 dj(r]), expressed by Wigner rotation matrices and normal- 
ized to dj(i] = 0) = 1. An effective GPD model at given input scale Qq is 
provided by taking into account three partial waves, 

F^t) = d j ( V )f j + \t)+r ] %. 2 (v)f j - \t)+ V i d j . i {r,)fr\t), (11) 

valid for integral j > 4. Such a model allows us to control the size of the 
GPD on the cross-over line and its Q 2 -evolution, see right panel in Fig. 3. 



August 18, 2010 0:57 



WSPC - Proceedings Trim Size: 9in x 6in JLABproc 



3. Extracting CFFs and GPDs from DVCS measurements 

Certainly, the access to GPDs from experimental data requires some soft- 
ware tools. The variety of both observables and models suggests setting up 
a flexible architecture, which would allow easy implementation of new mod- 
els, processes, data sets and fitting strategies. We wrote software prototypes 
based on object-oriented programming paradigm (using Python) and, al- 
ternatively, on functional programming (using Mathematica). In both cases 
we group scattering processes, theoretical frameworks, models, and exper- 
imental data in classes, which also serve as databases. The cross sections 
are implemented as functions of electromagnetic form factors and CFFs, 
where models for them are set up separately. Experimental data are stored 
in ASCII files, similar to commonly used ones but more standardized. They 
contain all information needed for the evaluation of observables. After spec- 
ifying data files, models, and conventions, our prototype software provides 
theory predictions for data points, depending on model parameters, which 
can be controlled by a fitting routine. 

In a first global fit 31 to photon electroproduction off unpolarized proton 
we took sea quark and gluon GPD models with two SO(3) partial waves 
at small x, reparameterized the outcome from HI and ZEUS DVCS fits at 
Q 2 = 2 GeV 2 , and employed it in fits of fixed target data within the scaling 
hypothesis. Thereby, we used the dispersion relation (4,5), where 

-«(t) 

l + x\l + x) \l + x] V" l + xM 2 / 

specifies a valence-like GPD on the cross-over line. Here, the skewness ratio 
r = lim^o H(x, x)/H(x, 0), a(t) = 0.43+0.85 t/GeV 2 , b controls the x ->• 1 
limit, and M the residual i-dependence, which we set to M = 0.8 GeV. The 
subtraction constant is normalized by d and Md controls the i-dependence: 

m=d(i-^ l y 2 . d3) 

We also included the parameter-free pion-pole model for the E GPD 32 and 
parameterized the H GPD rather analogously to Eq. (12) with b = 3/2. 

For the fixed target fits we chose two data sets. The first contains twist- 
two dominated (preliminary) beam spin asymmetry A^g and beam charge 
asymmetry A^ c coefficients from HERMES 12 ' 33 and 12 beam spin asym- 
metry coefficients A^g , which we obtained by Fourier transform of selected 
CLAS 35 data with small —t. The second data set includes also Hall A mea- 
surements 36 for four different t values. In light of the discussion 37 of Hall 
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Fig. 1. Experimental vs. model results, 31 pinned down by fits without (squares) and 
with (circles) HALL A data, for fixed target kinematics: Agg (1-18), A^q (19-36), Ag^, 
(37-54) from Ref. 34; A^g (55-66) and S^'™ (67-70) are derived from Refs. 35 and 36. 

A data, we projected on the first harmonic of a normalized beam spin sum 

eW'™= r^cosw)- — r d W ~ — fv ~ , (Mi 

BS i ^ dx Bj dtdQ 2 d(f> J J dx Bi dtdQ 2 d(j)' v ' 

where dw oc Vi{4>)V2{4 l )d(t) includes the Bethe- Heftier propagators, and we 
also neglected then the helicity dependent cross sections (beam spin differ- 
ences) . We confirm that formed beam spin asymmetries are compatible with 
CLAS ones 35 and we spell out that the second harmonics in HALL A data, 
i.e., effective twist-three contributions, are tiny or hard to separate from 
noise. Such contributions are small d in HERMES kinematics, too, where 
the constant A^,, appearing at twist-three level, is a twist- two dominated 
quantity that, as expected, 11 turns out to be correlated with A^. 

To relate the CFFs with the observables we employed the BKM for- 
mulas 11 within the 'hot-fix' convention. 38 First we excluded the Hall A 
data and set H to zero. A least squares fit (x 2 /d.o.f. « f) provides the 
parameters 

r val = 0.95, fo val = 0.45, d = -0.24, M d = 0.5GeV, (15) 

from one local minimum, cf. description of recent data in Fig. 1 (squares). 
This should not be considered a unique solution; however, it is compatible 
with the expectation that the skewness effect at small x should be small, 
i.e., r ~ 1, and that, according to counting rules, 39 b should be smaller 
than the corresponding j3 value of a PDF. The smallness of 6 val indicates a 
rather strong enhancement effect in the resonance region. In the second fit 
we included Hall A data, which remains challenging due to the steepness of 
data points (-0.33 GeV 2 < t < -0.f7GcV 2 ) — last four points in Fig. f. 



d Except for 3x2 beam spin asymmetry data points at large —t, xb, and Q 2 . 
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Fig. 2. ZSmH/ir obtained from different strategies: our model fits 31 (dashed (solid) 
curve excludes (includes) Hall A data), seven- fold CFF fit 42 with boundary conditions 
(squares), H, H CFF fit 41 (diamonds), smeared conformal partial wave model fit 43 within 
H GPD (circles). Circles (diamonds) are slightly shifted to the left (right) hand side. 
The triangles result from our neural network fit, cf. Fig. 3 (left). 



Again we took some local minimum with x 2 /d.o.f. ~ 1, giving 



„val 



1.11 



ival 



2.4, 



d 



-6.0, M 



sub 



1.5 GeV. 



(16) 



One observes a slight increase of r and a larger value of b , i.e., the 
enhancement of the GPD H in the resonance region diminishes. In agree- 
ment with a chiral quark soliton model estimate, 40 the subtraction constant 
remains negative and is now sizable; however, at present a positive sign 
cannot be excluded. A closer look reveals that our first fit fails to describe 
Hall A beam spin sums and underestimates the beam spin differences by 
about 50%, while our second one still underestimates all the cross sections 
by about 25%. For the latter we find a rather large remainder, effectively 
parameterized by H , which is roughly five times bigger than expected. Lon- 
gitudinally polarized target data provide a handle on H, 11 where CFF fits 41 
in JLAB kinematics provide at the means a two to three times bigger H 
contribution compared to our expectations (rg ~ 1, 6^ ~ 2). These findings 
are one to two standard deviations away from our big H ad hoc scenario. 

So far we did not study model uncertainties or experimental error 
propagation, since both tasks might be rather intricate. To illuminate 
this, we compare in Fig. 2 our outcomes for SmHfxB, t)/n versus xb at 
t = —0.28 GeV 2 (left) and for Hall A kinematics xb=0.36 versus —t (right) 
with results that do provide error estimates. The squares arise from con- 
strained least squares fits 42 at given kinematic means of HERMES and 
JLAB measurements on unpolarized proton, where the imaginary and real 
parts of twist-two CFFs are taken as parameters. Note that 3m£ and the 
other remaining eight CFFs are set to zero, however, all available observ- 
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Fig. 3. Neural network extraction of 5ReW(x B j, t)/n from BCA 34 and BSA 35 data (left). 
Ratio r = H (x, x)/H(x, 0) of LO sea quark and gluon GPDs for small x versus Q 2 (right). 

ables, even those which are dominated by these CFFs, have been employed. 
This might increase 'statistics', however, yields also a growth of systematic 
uncertainties. The huge size of the errors mainly shows to which accuracy 
H can be extracted from unpolarized proton data alone. 11 A pure H GPD 
model fit 43 (circles) to JLAB data provides much smaller errors, arising 
from error propagation and some estimated model uncertainties. Both of 
our curves are compatible 6 with the findings 42 and the H GPD model anal- 
ysis 43 of CLAS data. However, for Hall A kinematics the deviation of the 
two predictions that are based on H dominance hypothesis, see dashed 
curve and circles in the right panel, are obvious and are explained by our 
underestimation of cross section normalization of about 50%. Moreover, the 
quality of fit 43 % 2 /d.o.f. ~ 1.7, might provide another indication that CLAS 
and Hall A data are not compatible within this hypothesis, see, e.g., the 
two rightmost circles in the left panel for CLAS (xb = 0.34, t = — 0.3GeV 2 , 
Q 2 = 2.3GeV 2 ) and Hall A {x B = 0.36, t = -0.28GeV 2 , Q 2 = 2.3GeV 2 ). 
The pure T-L and H CFF fit 41 (diamonds), including longitudinal polarized 
target data, is within error bars inconsistent with the H dominated sce- 
nario (circles), however, (accidentally) reproduces our dashed curve. All 
of these exemplifies that within (strong) assumptions and the present set 
of measurements the propagated experimental error cannot be taken as an 
estimate of GPD uncertainties. An error estimation in model fits might be 
based on twist-two sector projection technique, 11 boundaries for the super- 
ficial model degrees of freedom, and error propagation. 

Let us give a short outlook. Neural networks may be an ideal tool to ex- 



e Note that in all fits the unpolarized HALL A cross section at —t = 0.33 GeV 2 is not 
well described, see fourth to the last data point in Fig. 1 and the rightmost square in 
the right panel in Fig. 2, which results from a x 2 /d.o.f. rs 2.3 fit. 



August 18, 2010 0:57 



WSPC - Proceedings Trim Size: 9in x 6in JLABproc 



10 

tract CFFs or GPDs. We present in Fig. 3 a first example in which, within 
H- dominance hypothesis, H is extracted by training 50 feed-forward nets 
with two hidden layers to HERMES BCA 34 and CLAS BSA 35 data. Hence, 
only the experimental errors were propagated, which in absence of a model 
hypothesis get large for the t — > extrapolation. Furthermore, we note 
that our Fortran code, used so far for small x, is now combined with the 
dispersion relation model for valence quarks, which gives us the possibility 
to provide predictions for the photon electroproduction cross section over 
a wide kinematical range. Our new fitting results are compatible with the 
presented ones, where the GPDs are modeled with three effective SO (3) par- 
tial waves (11). Thus, to LO accuracy the skewness ratio for gluons is now 
r ~ 0.5 at small x and rather stable under evolution, cf. Fig. 3 and Fig. 4 
of Ref. 31. Finally, we would like to add that the perturbativc description 
of hard exclusive meson electroproduction data is under investigation. 46 

In conclusion, a phenomenological access to the proton spin sum rule can 
be only reached within an understanding of GPD models, which can be also 
formulated in terms of an effective nucleon LCWF. A whole framework is 
available to reveal GPDs and to access the nucleon wave function, which 
should be considered as the primary task. To do so, reliable data, some 
mathematical understanding, and appropriate software tools are required. 
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'Note that this rules out the small- a: claim, 44 based on incomplete considerations at LO, 
which would roughly give for quarks and gluons r ~ 1.6 and r ~ 1.1, respectively. The 
reader might find further details in Ref. 45. 
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